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TECHNICAL MEMORANDUM X-64768 


SPACE SHUTTLE RENDEZVOUS, RADIATION 
AND REENTRY ANALYSIS CODE 

SECTION I. INTRODUCTION 


The austerity of the total space budget and the need to continue meaning- 
ful space-oriented scientific research demanded the development of a reusable 
Space Shuttle System. Comprehensive and detailed mission analysis for such 
a system is required to provide the kind of vehicle design and mission inte- 
gration data needed during Phase A and B studies and beyond, to ensure the 
development of a Space Shuttle System that performs as intended. 

The purpose of this effort was to develop a comprehensive, unified and 
versatile mission design and analysis tool which could be economically used to 
provide the kind of data mentioned above for preliminary investigations of the 
Space Shuttle System. This document may also serve as a ready reference for 
near earth orbital analysis and various parametric studies. 

The underlying philosophy of approach in developing this computer pro- 
gram was to produce a self-contained, multipurpose, multioutput package, 
placing emphasis on accessibility and usability with a minimum of input prepa- 
rations, while constraining accuracy within acceptable tolerances, thereby 
economizing with regard to computer run time. 

The Space Shuttle Rendezvous, Radiation and Reentry Analysis Code 
has been developed, checked out and used to provide some preliminary data 
for proposed earth resources Sortie Lab experiments. This document is not 
intended as a user's manual but rather to describe the mathematical model 
used and to give insight on the applicability of this Mission Analysis Code to a 
certain class of problems involving the Space Shuttle System. A complete 
program listing is available upon request. 

The Space Shuttle reentry analysis includes giving instantaneous reentry 
and landing site acquisition constraint parameters, all as a function of a 
specific range of atmospheric reentry angles and the reentry target altitude. 
Thus, suitable reentry conditions are achieved for the Orbiter at the beginning 
of reentry proper; i. e. , at the start of its plunge through the dense atmosphere. 



Section II gives a complete description of the Space Shuttle Mission 
Analysis code, including the mathematical modeling of the systems geometry 
and equations. Also, specific input-output options are detailed. 

Section m presents data generated using the Space Shuttle Mission 
Analysis Code to perform rendezvous radiation and reentry analysis, simulating 
five examples of Space Shuttle missions. It is divided into three major parts. 
Part A gives the results of parametric two-dimensional study on single impulse 
versus double impulse deorbiting for the Space Shuttle with constraints in the 
range of allowable atmospheric reentry angles and velocities at the reentry 
target altitude due to vehicle heating considerations. Also, the relative orbital 
manuevering system AV requirement is parameterized for deorbiting from 
varying altitudes and achieving reentry angles over the indicated range. Part 
B illustrates the results of performing three-dimensional mission analysis for 
five typical space Shuttle missions including a version of the proposed Apollo/ 
Soyuz Rendezvous and Docking Test Mission. Part C shows the results of the 
space radiation analyses of magnetically trapped protons and electrons encoun- 
tered during specific missions simulated in Part B. 


SECTION II. PROGRAM DESCRIPTION 
(SYSTEMS GEOMETRY AND EQUATIONS) 


Methods of orbital analysis and trajectory generation are many and 
varied. For preliminary analysis it is usually sufficient to perform basically 
a two-body Keplerian analysis with minimum effects, due to perturbing forces 
of some prescribed nature. There is a direct transformation between the six- 
dimensional space of Cartesian coordinates X, Y, Z, 5c, Y, fc and the six 
classical orbital elements a, e, t, ft, co, and i. 

A basic approach to trajectory generation, then, is to define an initial 
set of orbital elements along with any pertinent or secular time rates of changes 
of the elements, thus iterating between the two frames of reference performing 
the desired mapping of position and velocity at any instant of time. 

This section illustrates the coordinate system geometry and defines 
the equations used in constructing the mathematical models employed by the 
Space Shuttle Rendezvous, Radiation and Reentry Analysis Code to perform 
orbital analysis and trajectory generation for the Space Shuttle System. 
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In some instances only general geometric and functional relationships 
are presented because the detailing of specific manipulations is only pertinent 
to the actual programming effort. 

A. Basic Two-Body Geometry and Dynamics 

1. Orbital Plane Geometry. By definition of the properties of an 
ellipse, the following basic relationships are immediately derived from Figure 
1, assuming the variables a, e and 0 are available inputs: 

r = a(l-e 2 )/l+e cos 0 (l) 

p = a( 1-e 2 ) (2) 

b = a(l-e 2 )^ (3) 

where, by classic definition 

e = eccentricity 

a = semimajor axis 

b = semiminor axis 

p = semilatus rectum 

0 = true anomaly 

r = radius vector from focus 

From equation (l), with 0=0° and 180° respectively, we obtain 

r = p/(l+e) = a(l-e) (4) 

P 

r = p/( 1-e) = a( 1+e) (5) 

3 , 

which gives functional values for the radius of perigee and the radius of apogee, 
respectively. 
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Again from Figure 1 expressing the rectangular coordinates (X q , Y ) 

of the orbital plane as a function of the eccentric anomaly E and the true anom- 
aly 0, gives rise to the following relationships between E and 0: 


cos E 


e + cos 0 _ a - r 

1 + e cos 0 ae 


and 


sin E 


( 1-e 2 ) ^ sin 0 
1 + e cos 0 


( 6 ) 


(7) 


The focal radius equation thus becomes 


r = a( 1-e cos E) 


( 8 ) 


The total area A of the ellipse integrated over E gives rise to an angular 
parameter mean anomaly M 


M = E - e sin E 


(9) 


which is defined as a central angle compared to a circle having the same total 
area as the ellipse. 

2. Central Force Field Dynamics. The areal velocity or angular 
momentum vectors with the magnitude expressed in polarform reduces to 

i_ 

1 A l = y<i Ir* v| = ^ r 2 © = — = constant (10) 


This property permits the differential equation of motion to be easily 
derived, thus establishing the following parametric relationships 


P = 


h 2 


( 11 ) 
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from which 


h = [ a ( 1 -e 2 ) ^ ^ (12) 

where fj. = Gm, G being the universal gravitational constant and m the mass of 
the central body. 

Comparing the properties of equation ( 10) to the total integrated area 
of an ellipse and using previously defined quantities, we obtain the general 
time of flight equation 


t 


— (E-e sin E) 


(13) 


Thus, the total time for one revolution of an orbit is 


T = 



(14) 


The mean motion n defines the average angular rate of a body in orbit 
and is given from equation ( 14) to be 


n 



(15) 


where the units of n are radius per unit time. 

The mean anomaly M can now be defined also as a function of mean 

motion 


M = n(t-T) 


( 16 ) 


where t is the time of perigee passage. 


6 



Comparison of equation (9) and (16) yields 


M = n(t-r) = E-e sin E 


( 17 ) 


Thus, M is readily obtained when (t-r) is given and a good approxi- 
mation of E can be found when M is known from the following series 
expansion 


E 


M + e sin M + 


e 2 sin 2 M 
2 




(18) 


A general expression for the instantaneous linear speed of the orbiting 
body is derived from equation ( 10) to be 


V 


h 

r siny 


(19) 


where y is the angle between r and v and is called the flight path angle. 

3. Two Body Energy Relationships. A very useful computational energy 
parameter will be introduced and defined as 


A = 


r v 2 


( 20 ) 


where A is twice the ratio of kinetic to potential energy of the orbiting body. 
The semimajor axis may be defined now as 


( 21 ) 

2 - A 


and the eccentricity becomes 
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e = n/ 1 - A ( 2 - A) sin 2 y 


( 22 ) 


also the semilatus rectum can be written as 


h 2 r 2 v 2 sin 2 y , . 2 

p = — = = r A sim y 


(23) 


Thus, it is seen that the orbit-shaping parameters may be determined 
as a function of A andy. Using the above relationships, still another focal 
radius equation may be written: 


r 


r A sin 2 y 
o o o 


1 + vl - A Q ( 2 - A q ) sin 2 y Q cos © 0 


(24) 


Equation (24) is a very practical tool for computing the range of 
ballistic trajectories by treating them as fictitious orbits about the earth. 

Applying certain initial conditions along the major axis of the orbit as 
shown in Figure 2, some useful computational relationships become evident. 
Perigee injection equation (24) becomes 


r A 

— r P 

r l _ \J 1 _ A ( 2 - A ) 
P p p 


Solving equation (25) for A , we obtain as one solution 


2 r (r - r ) 

A = a P a 

P p 2 p 2 

P a 


(25) 


( 26 ) 


which defines A as a function of desired r . 

P a 


F or the circular orbit 
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Figure 2. Coplanar rendezvous and reentry geometry. 


A = — = 1 

P c r p 


( 27 ) 


Proceeding in a similar manner for apogee injection we have 


A = 


2 r (r - r ) 
P P a 

r 2 - r 2 
P a 


(28) 


defining also A as a function of the desired r . Again for the circular orbit 
a p 

case 


A = -? = 1 
a r 
c a 


(29) 


the ratio of A^ to A^ reveals a useful relationship for tying down an orbit. 
From equation (26) and (28) we obtain 


r 

A = — A 
p r a 


(30) 


and 


A = — A 
a r p 


(31) 


This gives the capability to determine the energy parameter value at 

the opposite end of the major axis of a specific orbit when the value for either 

end is known. Also the energy parameter values A and A are directly re- 

p a 

lated to the satellite velocities at each point, i. e. , at perigee V and at apogee 

V . P 

a 
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4. Time of Flight Computation Along Reentry Path. The time of flight 
parameter from deorbit retro fire to some specified reentry target point as 
shown in Figure 2 will vary depending on the desired atmospheric reentry 
angle y* and the altitude of the deorbit manuever. The specific problem, then, 
is to define the time it takes to go from point F to point G, or from point J to 
point K of Figure 2 after retro fire. 

Using the equations of Section I, we can express the time it takes to go 
from point H to point G of the deorbit conic shown in Figure 2 as 


t 


D 


P Ml 
D D 


2 7T 


(32) 


where P and is the period and mean anomaly of the deorbit conic respec- 
tively. Therefore, the time of flight from point F to point G may be expressed 
as 


T = — - 
D 2 


P^ M 

_D D 

2 7T 


(33) 


which reduces to 


t d 



(tt _ M d ) 


(34) 


where 


is computed from equations (6) and (17) 


B. Orientation of the Orbital Plane 
and Trajectory Generation 

The orientation parameters are designated S2, longitude of the ascending 
node; co, the angle or argument of perigee and i, the inclination of the orbit. 
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Figure 3 shows the particular frame of reference used in the Space 
Shuttle Mission Analysis Code to describe the relative position of the 
spacecraft. 

The reference longitude is the Greenwich meridian point G or the 

o 

Greenwich hour angle at the time of launch or insertion. The time-varying 
ascending node angle £2* is indicated as a function of rotational rate of the 
earth, and the nodal regression rate £2 of the particular orbit. 

£2*(t) = £2* + (SJ e + ft) t (35) 


where £2* is the instantaneous value at t = 0; i.e. , at launch or insertion and 
o o 

t is the elapsed time since t . 

o 

The position and velocity components of the system as a function of time 
are thus described by performing the indicated axes rotations of Figure 3 which 
is a transformation from orbital coordinates to geocentric Cartesian coordinates. 

Iterative or successive use of these transformation equations with the 
time-varying orbital elements serve to generate desired Space Shuttle 
trajectories. 

The three indicated rotations will be performed in the following order: 

(l) a positive rotation through angle £2* about the z-axis (2) a positive 
rotation through angle i about the new x-axis, i.e., X"', and (3) a positive 
rotation through angle a) about the new z-axis, i.e., Z". 

The matrix operation that completes the transformation from geocentric 
to orbital coordinates is then 


0 


1 cos c 0 

sin a? 

1 

Y 

__ 

-sin o) 

cos c 0 

0 

0 




1 

1 Z J 


i 0 

0 

J 

\ °l 

\ 


1 


A 0 0 U os £2* sin £2* °\ / x g\ 


0 cos i sin i 


-sin £2* cos £2* 0 


yO -sini cos y y 0 


0 1/ \Z 


G 


(36) 


/ w 
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Figure 3. Geocentric reference coordinate system. 




Performing the above indicated matrix operation yields 



where the elements of matrix equation (37) are defined to be: 

A tl = cos oo cos £2* - sin c o cos i sin £2* 

A 12 = cos co sin £2* + sin c o cos i cos £2* 

A l3 = sin co sin i 

A 2 i = - sin co cos £2* - cos co cos i sin £2* 

A 22 = - sin co sin £2* + cos co cos i cos £2* 

A 23 = cos co sin i 
A 3 i = sin i sin £2* 

A 32 = sin i cos £2* 

A 33 = cos i 

To solve for geocentric coordinates we need the inverse form of our 
rotation matrix which for orthogonal coordinate systems is just the transpose 
of the rotation matrix. 

Substituting also the derived relationships between the orbital plane 
coordinates (X , Y ) and the eccentric anomaly E from Figure 1, we obtain 

for geocentric rectangular coordinates: 
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Differentiating each of the three component equations of matrix equation 
(38) with respect to time and solving for E by differentiating equation (13) with 
respect to time, we obtain the velocity components in geocentric coordinates as 
a function of eccentric anomaly E. 



Position and velocity components as a function of the true anomaly 9 
also fall out directly in a similar manner, yielding 



and for the velocity components 
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Since we have computed the position of the orbiting body in rectangular 
coordinates, positions in spherical coordinates (r, $, A) are obtained from the 
transformation shown in Figure 4: 


r = + Y 2 + Z 2 


sin $ 


Z 

r 


tan A 


Y 

X 


(42) 

(43) 


(44) 


with the inverse relation being 


X = r cos $ 

cos A 

(45) 

Y - r cos 

sin A 

(46) 

Z = r sin 


(47) 


C. Orbit Determination From Initial Launch Site and 
Orbit Insertion Conditions 

The methods of orbital analysis presented in Section II, Part B, 

assumed an initial set of orbital elements a , e , r . u , S2* and i . It mav 

o o o’ o o o J 

be, however, that the only initial conditions information available is ( l) launch 

site geocentric spherical coordinates (r , $ , A ) , (2) launch azimuth (3) 

L L L 

the desired (r, v , y) for orbit insertion and (4) the orbit insertion geocentric 
spherical coordinates (r^., 4*^, A^) . The problem is then to obtain an initial set 

of dimensional and orientation orbital elements from this set of given informa- 
tion. Certain indicated constraints placed on some of the insertion parameters 
serve the purpose of simplification of method; however, the removal of these 
constraints does not add unduly to the complexity of the problem’s solution. 
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Once the initial set of orbital elements is obtained, the methods of orbital 
analysis and trajectory generation presented in Section n, Part B, are then 
initiated. 

Referring to Figure 5, the spherical coordinates (r , $ , X ) are the 

L L L 

launch site coordinates. The unit vectors (i , i , k ) define the orbital axes 

o o o 

and the unit vectors (i , j , k ) define the geocentric coordinate axes. The 

G G G 

spacecraft is launched from the designated site with the velocity vector lying 
along the launch azimuth 0 with flight path angle y. All initial conditions are 
referred to the orbit insertion point; i.e. , at perigee. The flight path angle y 

TT 

will thus be — and the actual orbit insertion point will be a down-range angle 

0* from the launch site. It is further established by definition that T will lie 

o 

along r ; J will lie along the launch azimuth and k will lie along the launch 
-L O O 

site meridian pointing north. Planar flight is assumed and when the actual 
flight is not planar, adjustments are necessary in the initial defining orbital 
elements. The transformation matrix from geocentric to nodal launch coordi- 
nate is obtained from the indicated rotation of Figure 5: 

( 1) a positive rotation about k through angle X 

G 

(2) a negative rotation about the new j axis; i.e., j' through angle 

- $ , and 
L 

( 3 ) a positive rotation about i" through the angle (90-0) which joins 
the two coordinate systems to coincide. 

The resulting matrix operation will then be 



(48) 


which yields 
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Figure 5, Orbit determination geometry from initial launch conditions. 




j a ll 3 12 

a \ 
a 13l 




Y 

o 

= 

a 21 

a 22 

a 23 

y g 

(49) 




^ a 3l 

a 32 

a 33^ 



A 


where 

a n = cos $ cos A 

^ Xj Xj 

a 12 = cos# sin A 

Xj Xj 

a 13 = sin # L 

a 2 j = - sin /3 sin A - cos /3 sin $ cos A 

Xj Xj Xj 

a 22 = sin/3 cos A - cos /3 sin # sin A 

Xj Xj Xj 

a 23 = cos /3 cos # 

Xj 

a 31 = cos /3 sin A -sin/3 sin# cos A 

Xj L L 

a 32 = - cos (3 cos A - sin /3 sin# sin A 

Xj Xj Ij 

a 33 = sin/3 cos # T 

Xj 

Element a 33 of matrix equation (49) gives the direction cosine for the 
angle between k and k' , thus defining the inclination of the orbit: 

(jr 

Tt* "P = cos i = sin/3 cos# . (50) 

Cj Xj 

From Figure 5 the unit vector along the line of nodes in the equatorial 
plane may be defined as 

r ft* = cos ft* T+ sin ft* "J . (51) 


20 


This vector transforms into orbital axes by the use of rotation matrix 
equation (49) and we obtain 


k Q = ^31 cos SI* + a 32 sin 
which leads to the relationship 


Tan 


tan A. - tan fi sin 4> 

■Li 

l + tanj3 tan X sin 

ij Li 


(52) 


If the first nodal crossing is a descending node, the ascending node is 

simply 

n* = n* ± iso 0 (53) 

a D 

depending on whether is east or west of Greenwich. 

For a due east launch /3 = 90° and the nodal longitude becomes from 
equation (52) : 


Tan ft* = Tan (90° + A ) 

L 


(54) 


which reduces to 


fi* = 90° + A 

L 


(55) 


giving a simple relationship for descending nodal crossing as a function of 
launch site longtitude. 

From Figure 5 an angle co* is defined. The relationship between co* 

and co is 


CO = CO* - © 


at perigee, 0 is equal to zero and we have 


( 56 ) 
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and 


CO — CO* 


i • r o* : cos a;* 
o 


Using matrix equation (29) we can define an initial angle co 

cos co = cos $ T ( cos cos X + sin sinA,^ (57) 

I \ a I a I J v ' 

where A^ is the insertion point longitude. Using previously defined relation- 
ships to compute shaping elements a, e and r we are thus able to compute all 
necessary initial defining elements which leads to a starting state vector for 
the Space Shuttle Mission Analysis at a zero reference time. Defining the 
variation of these elements with time enables us to compute a ground trace of 
the satellite’s position at any time beyond t as depicted in Section II, Part B. 

D. Formulation and Application of Certain Perturbing 
Effects on Keplerian Two-Body Motion 

Since the Space Shuttle Mission Analysis Code is concerned with pre- 
liminary orbital analysis, only limited use will be made of general perturbation 
theory and application affecting the motion of a Keplerian orbit. 

However, we may define with these analyses such effects as those due 
to (l) the oblateness of the Earth (2) a low thrust force vector (3) atmospheric 
drag and (4) a potential produced by the presence of other bodies of significant 
masses. The problem thus becomes to define the time rates of change of the 
six orbital elements as a system of linear differential equation and the various 
forms of the equations depending on the nature and origin of the force. 

We may write the general differential equation of motion for a perturbed 
Keplerian orbit as 


™ ( 58 ) 
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where R is the particular disturbing function and V is the operator, 


V 0X 


9 T + 


- 2-7 + 

9y ] 


9_ k 
8 z k< 


If we designate a. as being any dimensional element, a, e, M and/3^ as 
being any orientation element, £2, gj , i, and use the fact that r = f(t, (3 ^) , 
we may write using the method of variation of arbitrary constants, 


?£v?Sr*.- # 


1 1 


and 


1, 2, 3 


vr 

V 8a. i “ 9)3. i 


i i 


(59) 


(60) 


The simultaneous solution of equations (59) and ( 6 0) produces the 
transformation between the six-dimensional Cartesian coordinate space and the 
defining orbital elements ar's and j3's as a function of time. 



[a , a;.] a. + [a. , j3.] /3. 
r 1 1 r i i 


8R 

0QT 

r 


and 


r = 1, 2, 3 



[/3 , a.] a. + 
r x i 







(61) 


(62) 


Where any expression of the form [or, /3] is defined as a Lagrangian 
bracket and is of the Jacobian form 


[ar,/3] = Yj 
xyz 


9(X,X) 

8(«.^) 
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Using the explicit relationships of equations (36) through (41) in Section 
II, Part B, to obtain the differential coefficients of the direction — cosines, 
we resolve all the indicated Lagrangian brackets of equations (61) and (62) in 
explicit forms. 

Making these substitutions we arrive at the six linear differential 
equations depicting the time rates of change of the defining orbital elements as 
a function of the particular nature of the disturbing function R. 


a = 


2_ 9R 
na 9M 


(63) 


e = 1 - | ( 1 _e 2 ) — (l e 2 )^ 2 \ — 

na 2 e j 8M ~ U ’ e } f 9 cu 


(64) 


M - n . ) 2* _ 2 (23) 

\na 2 e / ' 9 e na \9a/ 


M 


(65) 


a = 


1 9R 

1/ 9i 

na 2 (l-e 2 ) 2 sin i 


( 66 ) 


CO - 


-COS 1 


na 2 sin i ( l-e 2 ) ^ 


9R + / ( l-e 2 ) ^ \ 9R 
9 i 1 na 2 e / 9 e 


(67) 


di 

dt 


cos 1 


0R 

1/ | 9 co 
\ h 


l /. 


na 2 sin i (l-e 2 ) 2 / \ na 2 sin i (l-e 2 ) /z 


1 < 68 ) 


To complete the analysis as it pertains to a specific problem, we need 
only to define R in an explicit form and take the indicated partial derivatives. 

Now, having an initial set of orbital elements a,e,M,fl,u>.i 

o o o o o 

and knowing the time rates of change of these elements, we can generate a 
trajectory supplying the Cartesian coordinate components of position X, Y, Z, 
and velocity X, Y, Z using the iterative method as shown in Section II, Part B. 
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This iteratiye method applied in a step-wise fashion actually amounts 
to a numerical integration of equations (63) through (68). 

1. Trajectory Generation for Tangential Low-Thrust Circle-to-Circle 
Orbital Transfers. We can uniquely describe the disturbing function operator 
V R in equation ( 58) to be a low -thrust tangential force vector F, where as the 
solutions of equations (63) through (68) will contain the additional components 
of acceleration, thereby perturbing Keplerian two-body motion. It is easily 
shown that starting with an initial circular orbit and initiating a low -thrust 
force tangent to the orbit, the most profound change occurs in the semimajor 
axis a. Thus, we can concentrate this analysis, obtaining in explicit form 

— in equation (63) giving us the solution to a. To expedite this we define 


VR 


22t + 2£t + 22t;. r 

9x 9y 9x 


(69) 


where T? can be explicitly resolved to always lie along the velocity vector by 
the transformation shown in Figure 3 and the transformation equations (36) 
through (4l) . We have then 


9R 

9M 


V R • 


9r 

0M 



(70) 


The explicit form for is obtained by using equation (9) and the 
transformation matrix equation (37). Thus equation (63) finally becomes 

2 [r 2 e 2 sin 2 Q + a 2 (l-e 2 ) ] 
r 2 n(i_e 2 ) /2 ( 1 + e 2 + 2e cos ©) ^ 

where F^ is the magnitude of the tangential force. Using equation (lO) we have 

1/ 

da _ da d© _ na 2 (l-e 2 ) 2 da 
dF " d© dt" d© 


(71) 
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We can now transform equation (7l) and show a change in the semi- 
major axis for one revolution of the orbit due to thrusting in a tangential 
direction: 


da = F„ 


2[r 2 e 2 sin 2 8 + a 2 (l-e 2 ) I 
r 2 n ( l-e 2 ) ^ ( l+e z +2 e cos ©) ( na 2 (l-e 2 ) /2 j 


V, 


d0 (72) 


It can be demonstrated that a characteristic of equation (72) is that the 
effect on Aa of a change in the eccentricity is at most second-order or higher. 

Also for tangential low-thrust, eccentricity changes very slowly with 
time and the e equation derived from equation (64) shows that when e = o, 
e = o. Thus if we set e = o in equation (72) and perform the indicated integral 
we obtain a Aa for one revolution of the orbit 


Aa = 



(73) 


Dividing both sides of equation (73) by the period of one orbit we obtain 



At n 


(74) 


We, therefore, arrive at a value of a for an initial circular orbit and 
perturb a at selected steps along the orbit. At the end of one revolution, a new 
a is calculated and a is now perturbed at the new rate. 

Applying the trajectory generation methods described in Section II, 

Part B, we compute position and velocity in geocentric Cartesian coordinates 
while the altitude of the orbit is constantly changing. 

2. Atmospheric Drag and Low-Thrust Descent. Atmospheric drag 
produces the direct opposite effect on Keplerian motion as that produced by a 
tangential low thrust force along the velocity vector, since the drag force vector 
is directly opposite the direction of motion. 
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We thus describe the nature of the drag force vector 


c dV v2 

F m — — 

2m 


where 

= drag coefficient 

= cross-sectional area of orbiting body 
m = mass 
V = linear speed 
p = atmospheric density at a point 

- F may be substituted in equation (74) and applied in the previously 
T D 

described method obtaining a value for -Aa due to drag. 

A conceivable Space Shuttle/Tug mission would be to retrieve a payload 
from geosynchronous orbit via low -thrust solar electric propulsion. If we 
direct the thrust vector F directly opposite the velocity, thus creating a 

force of magnitude - | F | , we can lower our orbit to the desired altitude, 

using previously described techniques for circle-to-circle transfers. 

3. First Order Perturbing Effects Due to an Oblate Earth. We may 
proceed in a like manner by defining the components of the force vector F 
arising from the potential produced by an irregular shaped — nonhomogeneous 
sphere. The components are resolved as before, relative to the inertial system 
described in Figure 3. 

Using equations (66) and (67) , along with previously shown methods, 
we describe the effects of an oblate earth only on those elements that undergo 
secular perturbation, thus ignoring the periodic perturbations. 


( 75 ) 
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These elements and their first order form are: 


a 


- 3 / 2 J 2 n cos i 



( 76 ) 


and 


% J 2 n (2 - 5 / 2 sin 2 i) 



(77) 


where the earth-zonal harmonic term is J 2 = 1.0827 by 10 -3 and the earth’s 
radius is R = 6378. 160 Km. 

It is obvious from equations (76) and (77) that the nodal regression 
rate and the perigee regression and advancing rates are primarily functions of 
the inclination of the orbit for purely Keplerian considerations. 

The perturbation in the mean anomaly M is compensated for by using 
a slightly perturbed value of mean motion for Keplerian orbits. 


An 


“ % J 2 n 



( 3 ,/ 2 sin 2 i-1) 


(78) 


Thus, the oblate mean motion becomes 


n* = n + An 


(79) 


and the change in the Keplerian period due to the earth’s oblateness becomes 


27T 

n* 


(80) 


T* = 


For low thrust trajectory analysis these values are updated after each 
revolution. 


E. Cross-Range/Down-Range Computation 

To meet the Space Shuttle vehicle constraints with regards to landing 
site acquisition opportunities (LSAO), the Space Shuttle Mission Analysis 
Code computes the instantaneous cross-range and down-range at each point in 
the program that position is computed. 

Performing two of the indicated rotations of Figure 5; i. e. , a positive 

rotation through X and a negative rotation through -$ , a vector is defined 
L L 

pointing to the launch site. 


L = cos $ cos X i + cos $ sinX j + sin 4> k 
L L Li L L 


(81) 


Performing two of the indicated rotations of Figure 3; i.e. , a positive rotation 
through S2* and a positive rotation through i, a vector perpendicular to the 
orbital plane is defined. 


H = sin i sin O* i — sin i cos O* j + cos i k 


(82) 


The instantaneous angle from which cross-range is computed becomes 

P* = sin -1 ( L*H) (83) 

To compute the instantaneous down-range, a unit vector perpendicular 
to the plane containing L and H is computed 


- L x H 
N = 

II* hi 


( 84 ) 
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Now, the projection of L on the orbital plane is given by 


P = H x n 


An instantaneous vehicle position in the orbital plane R (X, Y, Z) is 
computed as described in Section II, Part B. 

The instantaneous down-range angle which gives the relative position 
of the spacecraft in the orbital plane to the projection of the launch site vector 
on the orbital plane is thus 


©* = cos -1 (R*P) 


( 86 ) 


F. Space Radiation Analysis 

Since electrons and protons are magnetically trapped about the earth, 
a representation of their distribution may be made based on the contours of the 
magnetic field lines and the magnetic field strength at a point in space. This 
was accomplished by using the B-L coordinate system developed by Carl 
Mcllwain [10], as depicted in Figure 6. The B coordinate denotes the magnetic 
field strength at a specified point in space and the L coordinate is the magnetic 
shell parameter that specifies the shell upon which the guiding center of the 
trapped particles is confined as it drifts around the earth. The L coordinate 
is approximately constant along a geomagnetic field line. 

Essentially, a three-dimensional space of latitude, longitude and altitude 
is transformed into a two-dimensional space of B and L which serves to more 
expediently construct a model radiation environment based on this method. 

For the proton environment, the omnidirectional integral flux spectrum 
may be represented by 


Ej - E 
E (B,L) 

J (>E, BL) = F(B,L)e (87) 
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where E t to E is the particular energy band of interest and E q is a spectrum- 

shaping parameter and a function of B and L. F (B, L) is the known intensity 
of the proton flux for a given energy at a specific point. Equation (87) then 
defines the integral spectrum on the integral number of particles greater than 
Ej in the spectrum. 

The units of J are protons/cm 2 -sec. 



To obtain the differential energy spectrum which gives the number of 
particles at specific energies in the spectrum, we differentiate equation (87) 
with respect to E and obtain 


- J' (E,B,L) 


F(B,L) 

E (BL) 
o 


Ej - E 

e o (b,l) 

e 


where the units of J' 


are protons/cm 2 -sec- meV. 


( 88 ) 
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A differential and integral spectrum may then be completed for any 
point in the space model. However, for a typical Space Shuttle Mission, time- 
averaged calculation may be of greater interest. 


The Space Shuttle Mission Analysis Code calculates this time-averaged 
data for magnetically trapped protons in the following manner: a proton flux 
F[B(t),L(t)] is computed at each orbital step of 5 deg, which may vary, of 
true anomaly for an entire 7-day mission. 


From equation (87) we may compute this time averaged, or mean value, 
of the proton flux > Ej as 


J (> E) 




E t - E 

E [B(t),L(t)J 

F [B(t),L(t)]e 


dt 


t 


o 


(89) 


When E = Ej equation (89) takes the form 


J (> Ej) 



t 


n 



F [B(t),L(t)l dt 


(90) 


which gives a time-weighted average of all the partials in the spectrum greater 
than the specified energy Ej. Now, instead of choosing a representative 
spectrum at a single point, a representative average spectrum is chosen based 
on the spacecraft’s encounter with the radiation environment during the entire 
mission. From equation (88) we also have the time averaged differential 

energy spectrum for a Space Shuttle Mission with a time duration of t to t . 

o n 


J (E) 


- XLf 

f F (B,L) V^TTLOT 

J E Q [B(t),L(t)] 6 

t 

o 


(91) 
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where the units of J are as previously indicated. 

A similar analyses with some variation is performed to model the 
magnetically trapped electron environment. 


G. Program Implementation - Start/Stop - 
Input/Output Control Options 

1. Generation of Starting State Vector. The Space Shuttle Mission 
Analysis Code requires as initial input a starting state vector specifying the 
initial orbit along with an associated ground elapsed time (g.e.t.) since lift- 
off or g. e. t. since initial orbit insertion. This state vector takes the form of 
a set of geocentric orbital elements or other initial condition information 
defined in Section II, Part C. If planar flight is assumed during ascent, the 
code will generate its own starting state vector based on the following 
information. 

a. Launch site geocentric latitude and longitude. 

b. Orbit insertion geocentric latitude and longitude for the initial 
Shuttle base line orbit. 

c. Launch or insertion azimuth or the inclination of the desired orbit. 

d. The desired final altitude for on-orbit operations and the conducting 
of experiments. 

It is apparent then, that the method of obtaining the starting state 
vector will depend upon the kind of information available and other character- 
istics of the ascent portion of the Space Shuttle flight. 

2. Specific Start/Stop Program Control Options . Once the starting 
state vector has been defined, the following starting and run time control 
options are available: 

a. The ability to start the analysis at any time into the mission; i. e. , 
it may be desirable to perform only deorbit and reentry analyses at the end of a 
seven-day mission. 

This means also that any segmented portions of a mission may be 
analyzed for any desired time increments. For example, such an analysis may 
be required to define all the possible landing site acquisition opportunities 
during the entire mission. 
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b. Designating a mission time cutoff or specifying the number of 
revolutions desired as a means of terminating the analysis of a specific 
mission. 


c. A continuous unsegmented analysis of mission parameters from 
insertion into the operational orbit through deorbiting and the achieving of 
specified reentry constraint conditions (defined later). 

3. Rendezvous Analysis Output. Space Shuttle Mission Analysis Code 
provides the following output parameters to expedite rendezvous analysis: 

a. A complete ground trace of the spacecraft's trajectory from orbit 
insertion through the deorbit maneuver to the reentry target point. 

Two advantageous features of the ground trace computational tech- 
niques are, ( l) there is no dependence on knowing the relative position of the 
vernal equinox with respect to Greenwich as a function of time and, (2) the 
code employs a tracking technique which constantly updates those orbital ele- 
ments that undergo secular perturbations due to first order oblate spherical 
terms and the rotational motion of the earth, thus eliminating the need for 
separate time-consuming integration techniques. 

b. The position and the velocity of the spacecraft in any desired time 
increments defined in geocentric rectangular coordinates (X, Y, Z) and 
spherical coordinates (y, 4>, A.) . 

c. Orbital Maneuvering System (OMS) AV requirements for coplanar 
transfer and circularization maneuvers including phasing and retrograde 
deorbiting. 

d. Time into the mission associated with each event using the initial 
insertion time or g.e.t. since lift-off as the zero time reference. 

e. Nodal regression rates and perigee procession or regression rates 
along with the instantaneous nodal crossing. 

Rendezvous Analysis involving two vehicles is accomplished by, (l) 
r unn ing the target vehicle's trajectory, (2) modifying it if necessary using 
phasing orbits to change the original ground track to the desired ground track, 
(3) running the pursuit vehicle's trajectory starting with an initial state vector 
based on an instantaneous set of orbital elements defined from the target 
vehicle's trajectory, and (4) adjusting the pursuit vehicle's lift-off time and 
phasing altitude to insure a desired initial-phasing angle that corresponds with 
the desired rendezvous and docking time and position. 
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There is no automatic optimization techniques inherent in the Space 
Shuttle Mission Analysis Code; however, an optimum, or most economical, 
rendezvous sequence is arrived at through the operator's analysis of various 
pursuit and target vehicles parameters and trajectories. 

4. Reentry Analysis Output. As was stated earlier, the Space Shuttle 
Mission Analysis considers all space a vacuum, therefore reentry analysis 
does not attempt to define the spacecraft's trajectory through the atmosphere 
proper. The code, however, does define certain reentry constraint conditions 
and shows the time into the mission when all these conditions are met for 
successful Landing Site Acquisitions Opportunities (LSAO). Specifically the 
output includes: 

a. A complete deorbit profile from which any desired reentry angle 
into the earth's atmosphere may be chosen (usually a range from 0.0 to -2.0 
because of vehicle heating constraints) . 

b. The position and velocity of the spacecraft at the reentry target 
altitude as a function of the particular angle of reentry and the altitude at which 
deorbit retro fire was initiated. A single reentry angle may be chosen to 
expedite a faster run time and less data print-out. 

c. Instantaneous vehicle cross-range and down-range distances to 
landing site during the entire mission, which is used to determine at what time 
deorbit maneuvers may occur to allow for a successful LSAO. 

d. A running abbreviated deorbit profile (having a suppressed range 

of parameters) for possible mission aborts, which outputs all reentry constraint 
information as a function of deorbiting (retrograde) at any point of any circular 
orbit during the entire mission. This unique feature of the Code is 
accomplished by the periodic or cyclic rotation of the deorbit conic to coincide 
with whatever position in orbit that the spacecraft happens to be. 

e. The deorbit trajectories in Cartesian (X, Y, Z) and spherical 
coordinates (y, 4?, X) defined as a function of the particular angle of reentry 
chosen. Included in this output is the time of flight from deorbit retro firing 
to the atmospheric reentry target altitude which is also a function of range of 
reentry angles considered. 

f. Orbital Maneuvering System AV requirements for retrograde de- 
orbiting, also a function of desired reentry angle and altitude of deorbit initia- 
tion. There is no restriction here to circular orbits. However, when deorbit- 
ing occurs from other than a circular orbit it occurs at the apogee point of the 
orbit. 
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Again there is no attempt at automatic optimization, but an output which 
includes a full range of mission and vehicular constraint data enables the user 
to easily define all acceptable conditions of reentry, including the optimum 
ones. 


5. Radiation Analysis Output . Space Shuttle Mission Analysis Code 
performs a comprehensive analysis of the space radiation environment of 
magnetically trapped electron and protons encountered on a specific Space 
Shuttle Mission. The analysis included in this document was performed with 
the latest available environmental data; however, this data is periodically up- 
dated as new radiation environment models are defined and distributed. The 
output of this portion of the code includes: 

a. The differential and integral energy spectra for magnetically 
trapped electrons and protons as a function of spatial coordinates. 

b. A time-averaged differential and integral energy spectrum for 
protons and electrons as a function of mission time; i. e. , after any integral 
number of orbits or for the entire mission. 

c. Total number of particles above a specified energy encountered 
on a particular mission. 

d. Flux intensities for electrons and protons at any defined point in 
space thus defining the particular configuration of the model environment as a 
function of particle energies. 

These time averaged energy spectra may then be used to calculate 
radiation doses, including crew skin doses, by transporting them through 
certain thicknesses of materials using available nuclear radiation transport 
and dose calculation codes. 

6. General . The Space Shuttle Mission Analysis computer program 
was coded in Fortran IV — Double Precision, and currently runs on the IBM 
7044 computer. Since computer run time is greatly affected by the print-out 
option, the program source decks are now in three parts. Rendezvous and 
Reentry Analysis comprise a single deck requiring approximately 12 000 core 
storage locations. Radiation Analysis (Electrons) and Radiation Analysis 
(Protons) comprise two separate program source decks, each requiring approxi- 
mately 26 000 core storage locations. 

Print-out suppression options are available when less information is 
desired, thereby greatly decreasing the amount of computer run-time required. 


SECTION III. APPLICATIONS (COMPUTATIONAL 
RESULTS, TASK 1, 2, AND 3) 


The primary purpose of this section is to demonstrate the appli- 
cations of the Space Shuttle Mission Code to perform Space Shuttle Mission 
Analyses using realistic and probable missions. Rendezvous, radiation and 
reentry analysis has been performed for four typical Space Shuttle Missions. 
The presented data also includes a parametric study on single impulse versus 
double impulse deorbiting AV requirements for the Space Shuttle Orbiter 
which may be considered general reference data. 


A. Task 1 - Two Dimensional Reentry Analysis 

The first study undertaken was to show the relative orbital maneuvering 
system AV requirement for single impulse versus double impulse deorbit 
maneuvers when the Space Shuttle orbiter is operating at relatively low alti- 
tudes and reentering the atmosphere over the range of angles from 0.0 to -2.0 
deg. The range of reentry angles is dictated by vehicle heating constraints. It 
was discovered that, for specific angles within this range, the AV requirement 
for single impulse deorbit was considerably higher than the AV requirement for 
double impulse deorbit when the spacecraft is operating at altitudes below 400 
Km. Thus, it is more economical to transfer to a higher orbit before deorbiting 
for reentry when the reentry angle is within the applicable range. 

Figures 7 through 10 show the relative AV requirements for single 
impulse versus double impulse deorbit maneuvers at various orbiter altitudes 
and achieving the indicated range of atmospheric reentry angles. 

Using the ’’arrowed" single impulse line as a reference, all double 
impulse readings above the line represent a AV saving for double impulse 
maneuvers, and all double impulse readings below the reference line represent 
a AV penalty paid for performing the double impulse maneuver and achieving 
the indicated range of atmospheric reentry angles. 

Figures 9 and 10 show that as the initial circular orbit increases in 
altitude, little or no AV savings occur by performing a double impulse 
maneuver over the range of considered atmospheric reentry angles. 

Figure 11 shows the relative velocities that the spacecraft would have at 
the reentry target altitude of 120.38 Km (65 N. Mi.) after specified transfer 
maneuvers and reentering the atmosphere within the designated range of reentry 
angles. 
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Figure 11. Reentry velocity after transfer and deboost 
from initial 185.2 km circular orbit. 


B. Task 2 - Three-Dimensional Rendezvous 
and Reentry Analysis 

Task 2 demonstrates the versatility and flexibility of the Space Shuttle 
Mission Analysis Code to perform three-dimensional mission analysis for 
Space Shuttle Missions covering a wide range of orbital parameters. 

The illustrated computational results will show a ground trace of the 
orbiter from insertion into the operational orbit through the deorbit maneuver 
including the deorbit trajectory and the arrival of the orbiter at a reentry tar- 
get altitude, meeting a predetermined set of reentry constraints. 

Variation in the set of reentry constraint conditions has little or no 
bearing on the analysis of a particular mission, since a wide-enough range of 
instantaneous reentry constraint data is an output of each mission analyzed. 

One has only to define this set for a particular mission. 

Unless otherwise noted, the particular set of reentry constraint condi- 
tions, which assures successful landing site acquisition opportunities (LSAO) 
for the simulated missions presented herein, are as follows: 

(1) 4475 ^ down-range < 6475 n. mi. 

(2) O.On. mi. < cross -range < 1100 n. mi. (east or west) . 

(3) Reentry angle = -1. 31 deg (Mission One) and -1.35 deg (Mission Two). 

(4) Reentry velocity = 25693 ft/sec (Mission One) and 25845 ft/sec 
(Mission Two) . 

(5) Maximum 7 day mission duration time. 

1. Example Space Shuttle Mission One (Earth Resources). The first 
example mission chosen is a projected earth resources technology Sortie Lab 
mission. Orbital parameters are chosen to provide maximum viewing time 
during daylight hours of the Chesapeake Bay region with an approximate center 
of 38.0 deg altitude and -76.0 deg longitude. 

The orbit altitude is 268.54 km (145 n. mi.) and is an approximate 
daily repeating orbit with an inclination of 89.73 deg. The launch and landing 
site is the Western Test Range (WTR) , Vandenberg Air Force Base, California. 

Table 1 shows the sequence of flight events for the mission from inser- 
tion into the initial 50-by-100 n. mi. orbit to the achieving of the final set of re- 
entry conditions as defined earlier. The specific LSAO for Mission One aborts 
and mission completion would meet the defined reentry constraints. Time- 
lining information is also a feature of Table 1. 
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TABLE 1. TYPICAL SPACE SHUTTLE MISSION ONE (ERTS) 
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TABLE 1. (Concluded) 
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Notes: (1) altitude: 268.54 km (145 n. mi.) (5) orbit selection: a 1-day repeating, near polar orbit insuring passes over 

(2) time of launch: for daylight passes over target the Chesapeake Bay region with approximate center of 

(3) inclination: 89.73 degrees 38. 0 deg latitude and -76.0 deg longitude 

( 4 ) launch azimuth: 182.0 degrees (6) launch site WTR 




Figures 12, 13 and 14 show the ground trace of the first 16 revolutions 
of the orbiter from circular orbit insertions at point A of Figure 12. Since 
the orbit is an approximate repeating one, the revolutions essentially represent 
the complete ground trace for the whole mission. 

Figure 15 illustrates at what points above the earth (geocentric latitude 
<f> and longitude A) retrograde deorbiting must occur for the indicated revolu- 
tion numbers in order to meet the previously defined reentry constraints for 
both mission abort and mission completion reentry conditions. 

Figure 16 shows the remaining landing site acquisition opportunities 
meeting defined reentry constraints. 

2. Example Space Shuttle Mission Two (Advanced HEAP Delivery). 

For Space Shuttle Mission Two, all orbit parameters including the launch and 
landing site will differ from those of Mission One. The Space Shuttle will be 
launched from the Eastern Test Range, Florida. The operational altitude will 
be 200 n. mi. and the inclination will be 28.5 deg, due-east launch. 

Table 2 gives the sequence of flight events and the Orbital Maneuvering 
System AV requirement for each maneuver performed to get on station and to 
deorbit. All LSAO are tabulated, meeting the previously defined mission and 
reentry constraint conditions for Mission Two. 

Figure 17 shows the ground trace of the Orbiter from insertion into the 
200 n. mi. circular orbit at point A for the first 16 revolutions. 

3. Example Space Shuttle Mission Three (COMM/NAV Geosynch 
Development) . Example Space Shuttle Mission Three requires a kickstage or 
tug to place a Communication/Navigation (COMM/NAV) satellite in a geosyn- 
chronous 35786. 1 km orbit with the option of achieving different positions 
along a longitudinal shift. The orbit parameters for the Space Shuttle will be 
essentially the same as for Mission Two; thus the LSAO will remain the same. 
From a final circular orbit of 200 n. mi. , the Space Shuttle will serve as a 
launch pad for final deployment of the COMM/NAV payload by a delta kickstage 
to the desired geosynchronous position. 

Table 3 gives the sequence of flight events along with the relative time 
increments for kickstage firing to achieve the desired hovering point in geo- 
synchronous orbit. 
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Figure 12. Typical Space Shuttle Mission One (ERTSL) revolutions 1-8 
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Figure 13. Typical Space Shuttle Mission One (ERTSL) revolutions 9-13 
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LANDING SITE ACQUISITION OPPORTUNITIES (LSAO) MEETING THE FOLLOWING REENTRY CONSTRAINTS: 

(1) 4475 nmi < DOWNRANGE < 6475 nmi (4) REENTRY ANGLE = -1.31° 1 394,947 5 

(2) 0 n mi < CROSSRANGE < 1100 nmi (EASTORWEST) (5) REENTRY VELOCITY = 25693 ft/secj AT C65nmi) 

(3) 7 DAY MISSION DURATION 
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TABLE 2. (Concluded) 
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Figure 17. Typical Space Shuttle Mission Two revolutions 1-15. 




TABLE 3. TYPICAL SPACE SHUTTLE MISSION THREE (PLACE COMM/NAV 
SATELLITE INTO GEOSYNCHRONOUS ORBIT - 28.5 DEG INCLINATION) 


SEQUENCE OF FLIGHT EVENTS 


Event 

Time 8 

of 

Initiation 

(Sec) 

A Time C 
to 

Next Event 
(Sec) 

Propulsion 

System 

Event 

AV, 

m/a 

(fps) 

Resultant 
ha hp 
km 

(n ml) 

Latitude 

(deg) 

Position 

Longitude 

(deg) 

Altitude 

km 

(n ml) 

Insertion into initial 
orbit of 

50 by 100 n. ml. 

0,0 

2617.9 

NA 

NA 

185.2 by 92,6 
(100 by 50) 

27,57 

-6G.71 

92.6 

(so) 

First 

Impulse 

burn 

2617.0 

2643.2 

OMS 

27,58 
( 90, 5) 

185.2 by 185.2 
(100 by 100) 

-27,74 

102.47 

185.2 

(100) 

Second 

impulsive 
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5261.1 


OMS 

54.02 

(177.2) 

370.4 by 185. 2 
(200 by 100) 



185.2 

(100) 

Third 

Impulsive 

burn 

7 963.2 

2758. 6 C 
8 275.9 
13 793.2 
19310.5 

OMS 

53.64 

(176.01) 

370.4 by 370.4 
(200 by 200) 

mi 

80.59 

370.4 

(200) 

Option no. 1° fire 
kick-stage achieve 
position no, l-(Fig,18' 

10721.8 

43 082.04 

kick-stage 

or 

tug 

£405.7 
(7 892.9) 

35 786.0 by 370.4 
(19 322 by 200) 


-110.73 

370.4 

(200) 

Option no. 1, fire 
kick-stage to circu- 
larize in geosync. 
pos. 1 

53 803.8 

NA 

kick -stage 
or 
tug 

1459.5 

(4788.5) 

35 786.1 
by 

35 786.1 

-27.79 

-111.48 

35 786, 1 
(19 322.9) 

Option no. 2 C fire 
kick-stage achieve 
position no, 2- 
(Flg. 18) 

16 239.1 


kick-stage 

or 

tug 

2405,7 
(7 892.9) 

35786.0 by 370,4 
(19 322 by 200) 

27.51 

-133,38 

370. 4 
(200) 

Option no, 2 fire 
kick-stage to circu- 
larize pos. no, 2 

50321.1 


kick-stage 

or 

tug 

1459.5 

(4788.5) 

35 786.1 
by 

35 786.1 

-27.79 

-134. 53 

35 876. 1 
(19 322.9) 

Option no, 3 C fire 
kick-stage achieve 
position no. 3- 
(Fig. 18) 

21756.4 


kick-stage 

or 

tug 

2405,7 
(7 892.9) 

35 786,0 by 370.4 
(19 322 by 200} 

27.41 

-156.04 

370,4 

(200) 

Option no. 3 fire 
kick-stage to circu- 
larize in geosync. 
pos. no. 3 

64 838.4 

NA 

kick-stage 

or 

tug 

1459.5 

(4788.5) 

35 786.1 

by 

35786.1 

-27.79 

-157,58 

35 786.1 
(19 322.9) 

Observe: 

Note: 

The time for kick-stage firing to achieve relative positions 1-15 shown in Tlgure 41 is sequenced exactly one revolution apart. 

All landing site acquisition opportunities for mission aborts and mission complete are the same as defined for typical Space 
Shuttle Mission Two. 

LSAO 108th rev. b 
final for 7 -day Shuttle 
Mission - deorblt 

601424.9 

1734.5 

OMS 

101.77 

(338.89) 

370,4 by 0,0 
( 200 by 0) 

-25,34 

JHHj 

370,4 

(200) 

Vehicle at 120.38 km- 
reentry angle -1.35 and 
all other constraints 
met 

603 159.4 

TBD 

NA 

NA 

370,4 by 0.0 
(200 by 0) 

-0.65 


120.38 

(65) 

Landing 

TBD 

NA 

NA 

NA 

NA 

28.5 


0.0 


a. time = elapsed seconds since insertion into 50 by 100 n. mi, orbit 

b. revolution = number of orbits since insertion into 200 by 200 n. mi, circular orbit 

c. A time * time Increment to use for firing of kick-stage or tug as a function of desired geosynchronous position 


Notes: 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


altitude: Shuttle - 370.4 km (200 n. ml.) - COMM/NAV - 357B6.1 km (1932.29) n. mi. 
time of launch: TBD 

Inclination: 28.5 deg (6) expendable or retrievable delta kick-stage or tug 

launch site - landing site: ETR (7) orbit selection: the 200 n. ml. Shuttle orbit may be lowered depending on 

launch azimuth: 90 deg final payload weight and propellant requirements. COMM/NAV is placed in a 

geosynchronous orbit with the option of selecting the hovering point over 
a 360 deg longitudinal shift as shown In Figure <8. 
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Figure 18 illustrates the relative positioning achieved by the COMM/NAV 
payload as a result of firing the kickstage at increments of one revolution 
(5 517.29 sec) of the Space Shuttle Orbiter operating in a 200 n. mi. circular 
orbit. 


As will be explained in detail later, Figure 18 also shows the relative 
time averaged magnetically trapped electron particle count greater than 0.5 
MeV at the different position in a geosynchronous orbit. 

4. Example Rendezvous and Payload Retrieval Mission Four. For 
Mission Four the Space Shuttle Mission Analysis Code performs a version of 
the proposed USSR -SOYUZ /USA -A polio Rendezvous and Docking Test Mission 
to demonstrate the Code’s capability to perform rendezvous analysis involving 
two vehicles launched at different times and possibly from different launch sites. 
The simulated mission as performed is similar to the joint project Technical 
Proposal [ 13] only in the fact that an effort is made to insure that major events 
occur at similar times and over similar points on the surface of the earth. 

Emphasis is also placed on minimizing the total Reaction Control 
System AV requirement for the Apollo pursuit vehicle in accomplishing the 
mission. 

A current set of proposed orbit parameters were used to generate 
starting state vectors for both the Apollo and Soyuz vehicles. The Apollo 
launch site is KSC, Florida (ETR) and the Soyuz launch site is Baikonur, 
Kazakhstan, USSR. Other orbit parameters for each vehicle are detailed in the 
following tables and figures. 

Table 4 shows a detailed sequence of flight events for the Apollo/Soyuz 
Rendezvous and Docking Test Mission including relative phase angles and A 
node angle for the two orbits. Many of the detailed onorbit operations are 
omitted, but again, time is allowed for the operations and is shown as "vehicle 
phasing". 

Figure 19 illustrates the six-impulse rendezvous maneuver sequence 
as used in the Space Shuttle Mission Analysis Code. 

Figure 20 shows the relative positions and revolutions of the Soyuz and 
Apollo vehicles at Apollo lift-off from KSC, Florida and the relative positions 
of the spacecrafts at the start of rendezvous phasing. 
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TABLE 4. USSR-SOYUZ /USA -APOLLO RENDEZVOUS TEST MISSION 


SEQUENCE OF FLIGHT EVENTS 


Event 

Time 8 

of 

Initiation 
Hr (Sec) 

A Time 
to 

Next Event 
Hr (Sec) 

Propulsion 

System 

Event 

AV, 

m/s 

(fps) 

Resultant 

hs/hp 

km 

(n mi) 

j Position 

Phase 

Angle 

(Deg) 

A Nude 
Angle 
(Deg) 

Latitude 

(Deg) 

Longitude 

(l^B) 

Altitude 

km 

( n mi) 

Soyuz lift-off 

0.0 

(0.0) 

0. 161 
(580.0)) 

NA 

NA 

NA 

45.38 

63.16 

0.0 

NA 

NA 

Soyuz Insertion 
101 by 123 n, mi. 

0. 161 
(580.0) 
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(2659.02) 

NA 

NA 
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NA 
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MB 



NA 
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113 by 113 n, mi. 

9.0 
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0.0086 2 
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tm 
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NA 
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NA 
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209.5 

19.03 

0.001 
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WEZBJt/k 

Apollo — phasing 
Soyuz — phasing 

10.50 
(37 802.83) 

1.46 
(5 342.73) 

NA 

NA 





17.51 

0.005 


BLUM 


227.8 

209.5 
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Soyuz — phasing 

11.98 
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1.48 
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W 1 ' I ^'4 , ? i ^ > - I>|1 
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15.99 

0.008 
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•MW ‘ , - 

13.46 
(48 4R8.29) 

1.48 
(5 342.73) 

NA 

NA 

209.5 by 209.5 

RBI 

-147.15 

209.5 


0.011 

IffltjEIlHQi 

49.58 


K2B3H 

Apollo — phasing 
Soyuz — phasing 

14.95 
(53 831.02) 

wmmm 

NA 

NA 


44.49 
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■H9 


0.015 


49.66 


227.8 

Apollo — phasing 
Soyuz — phasing 

16.43 

(59 173.75) 


NA 

NA 

209.5 by 209.5 

45.34 

171,83 

209.5 

11.44 

0.018 


49,73 

-172.47 
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Apollo — phasing 
Soyuz — phasing 

17.92 
(84 518.48) 

■ 

NA 

NA 
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46. 16 

151.43 

209.5 

9,92 

0.021 


49,80 

165.24 

227.8 

Apollo — phasing 
Soyuz — phasing 

19.40 

(69859.21) 

1.48 

(5342.73) 

NA 

NA 

209.5 by 209.5 

46.93 

131.08 

209.5 

8.40 

0.025 


49.87 


227.8 

Apollo — phasing 
Soyuz — phasing 

20.88 

(75201.94) 

1.48 

(5342.73) 

NA 

NA 


47.65 

110.81 

209.5 

6.88 

0.028 


49.94 


RS&R 

Apollo phasing 

Soyuz — phasing 

22.37 

(80544.67) 

1,48 
(5 342.73) 

NA 

NA 


48.32 


SKXJH 


0.031 


mam 


KEtlLR 

Apollo — phasing 
Soyuz — phasing 

23.85 

(85887.40) 


■91 

Ijfl^ 

209.0 by 209.0 

48.94 

70.45 

2 09.0 


0.035 


a 

□i 


Apollo — phasing 
Soyuz — phasing 

25.34 
(91 230.13) 

■B 



209.0 by 209.0 
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2.32 

0.038 


50.14 


227.0 

Apollo — phasing 
Soyuz — phasing 

26.82 
(96 572.8) 

0.051 

(186.35) 


Eli 

209.0 by 209.0 

50.00 

9991 

209.0 


0.041 



31.53 

227.0 

Apollo plane change 
Soyuz position 

26.87 
(96 7 59.15) 

■nsnfrav 

SPS 

6.11 
(20 . 04) 

209.0 by 209.0 

51.72 

49.40 

209.0 - 


0.0 


51.71 

50.61 

227.0 

Apollo (TPI) 
Soyuz position 

27.303 
(98 292.2) 


' 1 

HRS 


-10.80 

141.55 

209.0 


0.0 

(21.32) 



■E3HEM 

i 


Apollo (TPF) 
Soyuz position 

27.88 
(100 369.9) 

0.194 

(701.13) 

«g| 

6.5 

227.0 by 227.0 

9 


227.0 


0.0 

(21.32) 

227.0 by 227.0 

-19.98 

MELSm 

wmam 


28.07 
(101 071.04) 


RCS 

I&K3II 

227.0 by 227.0 

16.70 

Bi 


1 

0.0 

Apollo — Soyuz 
docking 

28.27 

(101772.1) 

TBD 

RCS 

2.2 

(7.21) 

227.0 by 227.0 

47.02 

-3.94 

227.0 

0.0 

0.0 


Notes: (l) Soyuz launch site: Baikonor, Kazakhstan, Initial orbit: 187. 1 by 227. 8 km, Inclination: 51.7227 deg 

(2) Apollo launch site: Kennedy Space Flight Center, Initial orbit: 149.9 by 183.3 km. Inclination: 51.7227 deg 

(3) Apollo Insertion condition: Launch azimuth: 39.38 deg; Latitude: 39.29 deg; Longitude: 65.08 deg; descending node: 155.259 deg 

(4) Apollo launch time: 7,3729 hrs after Soyuz lift-off 

(5) Apollo phasing orbit: 209.48 km circular (113 n. mi.) 

(6) Soyuz phasing orbit: 227.80 km circular (123 n, mi.) 
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APOLLO PHASING 
ORBIT 209.5 km 



SOYUZ CIRCULAR 
ORBIT 227.8 km 


1. INSERTION - 150 0 X 183 3 km 
2 IMPULSE BURN - 183.3 X 209.5 km 

3. IMPULSE BURN - CIRCULARIZATION 209.5 km 

4. PLANE CHANGE TO CORRECT A NODE ERROR 

5. TERMINAL PHASE INITIATION - 207 X 230 km 

6. TERMINAL PHASE FINAL - 227.8 X 227.8 km 


Figure 19. Orbit geometry of Apollo T s six 
impulse rendezvous maneuver sequence. 

Figure 21 illustrates the points above the earth where terminal phase 
initiation (TPl) and the final terminal phase (TPF) occur and the docking 
revolution which very closely approximates the time and position for these 
Soyuz/Apollo maneuvers as set forth in Reference 13. 
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USSR - SOYUZ / USA- APOLLO RENDEZVOUS TEST MISSION 
APOLLO PHASING ORBIT: 209.48 km CIRCULAR (113 n mi) 

SOYUZ PHASING ORBIT: 227.80 km CIRCULAR (123 n mi) 
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Figure 20, Apollo liftoff and begin Soyuz Apollo phasing revolutions* 



USSR - SOYUZ / USA-APOLLO RENDEZVOUS TEST MISSION 
APOLLO PHASING ORBIT: 209.48 km CIRCULAR (113 n mi) 
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-200 -160 -120 -80 -40 0 40 80 120 160 200 
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Figure 21. Apollo terminal phase and Soyuz Apollo docking revolution. 




C. Task 3 - Space Radiation Analysis 

The computation results of Task 3 will demonstrate the capability of the 
Space Shuttle Mission Analysis Code to perform space radiation environment 
analyses of magnetically trapped electrons and protons. These analyses will 
define the possible radiation hazards associated with a specific mission, thus 
providing useful data for consideration in mission design options. 

Elements of Task 3 included the integration of the codes' trajectory 
calculation techniques and space radiation environment models in the form of 
data decks provided by Dr. James I. Vette [9], et al. , containing the measured 
omnidirectional flux of trapped electrons and protons at all defined spatial 
points. 


Since there has been some expressed concern because of possible 
trapped radiation hazard associated with placing an advanced HEAO satellite in 
an initial orbit greater than 200 n. mi. , space radiation analysis for the 
example Space Shuttle Mission Two was performed for altitudes of 200 n. mi. 
and 500 n. mi. 

Time-averaged differential and integral spectrum data for this mission 
will be presented in subsequent figures along with comparative spectrum data 
associated with the two indicated altitudes, holding other orbit parameters for 
the mission constant. It is evident from the configurations of trapped electrons 
and protons in the South Atlantic anomaly that the time-averaged spectrum 
intensity is a function of the orbital inclination as well as altitude. 

Another aspect of defining the space radiation hazard involves the use 
of a solar electric propulsion stage as a final trig to transport and possibly 
retrieve payloads for certain Space Shuttle/Tug missions. 

The specific problem is then to define the number and energies of 
magnetically trapped protons and electrons encountered on a specific mission; 
converting these particles into 1 MeV -equivalent electrons and applying a 
tabulated solar cell damage factor to arrive at relative power loss for each 
mission. 

This can be done by flying a simulated low -thrust trajectory through 
the space radiation environment model and computing the particle fluxes at all 
spatial points along the way. The first simulation is flown without power loss, 
whereas a trajectory for comparison includes a power degradation model which 
is linear with respect to the particle accumulation rate. 
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By knowing the environment, the damage to the solar cells and, thus, 
the degradation of the available power may be assessed as a function of the 
particular mission. 

Obviously a loss of power lowers the available thrust of the stage and 
this fact must be taken into consideration when establishing time-lines for a 
particular Space Shuttle/Tug mission. 

When the power loss is known as a function of time, this enables a more 
accurate trajectory to be calculated, using previously defined methods. 

Experiments and theoretical analyses are currently underway to deter- 
mine the dose rates and damage factors on specific solar cell models after 
being bombarded by various energy levels of protons and electrons. 

Figure 22 shows the time-averaged proton flux encountered during 
typical Space Shuttle Mission Two (200 n. mi.) and the relative particle popu- 
lation at an altitude of 500 n. mi. for the same mission. 

Figure 23 shows the same information for time-averaged electron 
fluxes for Mission Two. 

Figure 24 gives the differential and integral energy spectra for protons; 
also the total number of particles encountered during the 7 days of Mission Two 
at 200 n. mi. 

Figure 25 shows the same information for protons, but at an altitude 
of 500 n. mi. , for Mission Two. 

Figure 26 shows the differential and integral energy spectra for electrons 
and the total number of particles of all energies encountered during the 7 days 
of Mission Two at 200 n. mi. 

Figure 27 shows the same information for electrons but at an altitude 
of 500 n. mi. for Mission Two. 

Figure 28 shows the proton isoflux contours (E > 50 MeV) in the South 
Atlantic anomaly at an altitude of 145 n. mi. 

Figure 29 shows the greatly expanded proton isoflux contours with 
higher intensities (E > 50 MeV) at an altitude of 500 n. mi. over the South 
Atlantic anomaly. 
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Figure 30 shows the electron isoflux contours for passes through the 
South Atlantic anomaly at an altitude of 145 n. mi. 

Figure 31 shows the projected effect of the magnetically trapped space — 
radiation environment on an unshielded solar electric low-thrust tug orbital 
transfer from a 20 000 km circular orbit to geocynchronous altitude at inclina- 
tions of 28. 5 and 0.0 deg. It should be noted, that for this particular 
mission, a relatively thin transparent glass shield would reduce the total 
power loss over the mission to less than 3 percent. However, if a mission is 
started at a significantly lower altitude, power loss due to solar cell damage 
will cause increased concern to mission analysts with regard to mission dura- 
tion and whether a particular mission can be flown. 

Figure 32 gives an indication of the kinds of problems associated with 
solar electric low -thrust orbital transfers if a mission is initiated at unaccept- 
ably low altitudes. With little or no shielding, the low energy proton environ- 
ment may cause accumulated power losses in excess of 50 percent over 
relatively short mission durations, which renders some missions impossible 
to be flown. However, adequate shielding of the solar cells can reduce the 
accumulated power loss to acceptable levels; i.e. , in the 10 to 20 percent range. 



500 X 500 n mi 



3.68 X 10 2 < AV. FLUX < 3.88 X 10 2 



TIME AVERAGED FLUX ENCOUNTERED FOR 28.5° INC 
AT ALTITUDES OF 200 AND 500 nmi. 





500 X 500 . nmi 



2.35 X 10 5 < A V. FLUX < 2.51 X 10 5 



TIME AVERAGED FLUX ENCOUNTERED FOR 28.5° INC. 
AT ALTITUDES OF 200 AND 500 nmi. 
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Figure 24. Typical Space Shuttle Mission Two 
time averaged energy spectra (protons). 
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PROTONS/cm -sec(mev) 
'applies to diff. spec, only 




ElECTRONS/cm - sec (mev) 
‘applies to di ff spec, only 



Figure 26. Typical Space Shuttle Mission Two 
time averaged energy spectra (electrons). 
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PARTICLE ENERGY, E (mev) 


Figure 27. Typical Space Shuttle Mission Two 
time averaged energy spectra (electrons). 
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Figure 28. Proton isoflux contours at an altitude of 268.54 km (145 n. mi.) South Atlantic radiation anomaly 
producing total proton flux > 50 mev (AP7 data) for typical Space Shuttle Mission One (ERTSL.) . 
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Figure 29. Proton isoflux contours at an altitude of 926 km (500 n. mi.) South Atlantic radiation anomaly 
producing total proton flux > 50 mev (AP7 data) for typical Space Shuttle Mission One (ERTSL) . 
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Figure 30. Electron isoflux contours at an altitude of 268.54 km (145 n. mi.) (1968 data, Vette) South 
Atlantic radiation anomaly producing partial electron flux >0.5 mev as seen by typical 

Space Shuttle Mission One (ERTSL) . 



SECTION IV. CONCLUSION 


It has been demonstrated that the Space Shuttle Rendezvous, Radiation 
and Reentry Analysis Code is a basic and versatile Space Shuttle Mission 
design and analysis tool which allows for extensive user interaction and 
flexibility through the utilization of a relatively small computer (IBM 7044). 

It is hoped that the foregoing presentations concerning the development and 
applications of the code will render insight to persons engaged in preliminary 
Space Shuttle mission designing, as to whether the illustrated features of the 
code are applicable to their specific studies and problems. The capability of 
the code is currently being expanded to include elliptical to circular low -thrust 
orbital transfers, a more sophisticated power degradation model and other 
effects on low -thrust trajectory analysis such as shadowing. 

Mission analysis data generated by this code have been compared 
favorably with data generated from other sources. Source decks, listings and 
other specific information concerning the use of the Code are available from 
the author upon request. 
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